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Hund's paradox and the collisional stabilization of chiral molecules 
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We identify the dominant collisional decoherence mechanism which serves to stabilize and super- 
select the configuration states of chiral molecules. A high-energy description of this effect is compared 
to the results of the exact molecular scattering problem, obtained by solving the coupled-channel 
equations. It allows to predict the experimental conditions for observing the collisional suppression 
of the tunneling dynamics between the left- and the right-handed configuration of D2S2 molecules. 
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Introduction. — An old problem in molecular quantum 
mechanics, first discussed by F. Hund [H, is how to ex- 
plain from first principles why molecules often appear 
as enantiomers , i.e., either in a left-handed configura- 
tion or as the right-handed mirror image. Given the 
parity-invariant molecular Hamiltonian, one might rather 
expect them in the ground state, corresponding to the 
symmetric superposition of these chiral states. Tradi- 
tionally, this is explained by the possibly very long tun- 
neling time from a left-handed configuration state \L) to 
a right-handed one \R). However, this does not solve the 
'paradox' since one still needs to understand the seeming 
failure of the superposition principle, prohibiting super- 
position states of the form |^) = +e"^|i?)) /y/2 
with s$ £ < 2 

While this super-selection phenomenon has been linked 
with fundamental parity violations a very natural 
explanation is offered by the concept of environmental 
decoherence @. What selects the enantiomer states ac- 
cording to this theory is the fact that the typical interac- 
tion with environmental degrees of freedom, such as the 
collision with a gas particle, can better distinguish the 
alternatives \L) and \R) than, say, between the molecu- 
lar eigenstates |i/>o) and The enantiomer states are 
then prevented both from tunneling between each other 
and from decaying into a mixture, if these environmen- 
tal interactions are sufficiently frequent compared to the 
tunneling period. This stabilization can be understood in 
analogy to the quantum Zeno effect if one views the envi- 
ronment as continuously monitoring the molecular state. 
Superposition states \ip^), in contrast, would get quickly 
decohered by their ensuing quantum correlation with the 
environment. 

This environmental distinction of specific molecular 
configurations is one of the paradigms in the field of de- 
coherence, discussed by many of the path-breaking works 
@, BS|- At the same time, not much is known about the 
concrete microscopic mechanism at work with realistic 
molecules, nor whether the transition from the tunneling 
regime to stabilization can be observed experimentally. 

In this letter, we use molecular scattering theory to 
identify the dominant microscopic mechanism responsi- 



ble for chiral stabilization due to a background gas. We 
show that this effect is determined by a parity-sensitive 
higher-order term in the dispersive interaction, which is 
usually disregarded because it does not affect the equi- 
librium properties of gases (due to the oricntational av- 
eraging involved). In spite of this, and in spite of the 
fact that this process can not be used to separate a 
racemic mixture into chiral components, it provides a 
surprisingly efficient channel for environmental decoher- 
ence. This is demonstrated numerically by solving the 
full-fledged coupled-channel problem for a simple chiral 
molecule, D2S2. Our numerical results motivate and con- 
firm a high-energy approximation, which allows to assess 
the decoherence effect at room temperature and to pre- 
dict when the collision-induced transition from tunneling 
to stabilized chiral states can be observed experimentally. 

Master equation for collisional stabilization. — The fol- 
lowing microscopic analysis is facilitated by the recent 
derivation of a master equation yielding the incoherent 
dynamics of the internal-rotational molecular state due 
to the collisions with a lighter, thermalized background 
gas @, [If • Crucially, this Markovian description incorpo- 
rates the interaction between molecule and gas particle 
in a non-perturbative fashion, by means of the multi- 
channel scattering amplitudes. 

We assume that the kinetic energy transferred in the 
collisions is sufficiently small compared to both the Born- 
Oppenheimer barrier separating the enantiomer states 
and to the excitation energies of the electronic and vi- 
brational internal states of the molecule and the gas par- 
ticle. This is valid in most cases of interest, and it im- 
plies that thermally induced transitions between \L) and 
\R) do not occur. It is then justified to take the inter- 
action operator to be diagonal in the enantiomer basis, 
V = V L (r) \L)(L\ + V R {r) \R)(R\, thus restricting the 
molecular configuration state to a two-dimensional sub- 
space. Here, f is the inter-particle position operator in 
the body-fixed molecular system; for simplicity we dis- 
regard a possible dependence of V on the orientation of 
the environmental gas particle. 

The master equation Q simplifies considerably under 
these assumptions. It is formulated in terms of the scat- 
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tcring amplitudes for the channels corresponding to the 
molecular eigenstates |?/>o) an d \ipi}- However, since the 
S-matrix for V does not couple subspaces of different 
chiral configurations, (L\S\R) = (R\S\L) = 0, one can 
express the proper scattering amplitudes as linear com- 
binations of the scattering amplitudes and f( R > as- 
sociated to the unitary scattering operators (L\S\L) and 
(R\S\R). The configuration dynamics then no longer de- 
pends on the scattering cross section. Rather, it is deter- 
mined by 
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which may be called a decoherence cross section. Similar 
to a proper partial cross section, this characteristic area 
depends on the relative velocity v and on the initial and 
final internal state of the molecule, labeled by the multi- 
indices ao and a. In the present case, these are the rota- 
tion states of an asymmetric top, a = (j, n%j, r), specified 
by the total and azimuthal quantum numbers j and rrij , 
and the pseudo quantum number r. Importantly, the de- 
coherence cross section (Q]) depends on the phase differ- 
ence of /( L ) and f( R \ It may thus be quite large even if 
the corresponding scattering cross sections are identical, 
|y(i)|2 _ |y(fl)|2_ rp-^g p resence f background gas 

also gives rise to a coherent modification of the tunneling 
dynamics, described by the characteristic area e aao (v) = 
J dndn Im[/a L 2 (vn,vn ) vn )]/ATr. 

In the master equation the areas rj aao (v) and e aao (v) 
are multiplied with the current density of the gas par- 
ticles to yield the decoherence rate 7 and the frequency 
shift u> x , respectively. Specifically, the decoherence rate 
is given by 

/>oo 

7 = n gas (v7]) = n gas ^2 w (ao) / dw (v) vr) aao (v) , 

a,ao 

(2) 

where n gas is the number density of the background gas 
and v (v) its velocity distribution. Here, taking the av- 
erage over the rotational state distribution w (ao) is an 
additional approximation, which is well allowed if the 
rotation frequency is much greater than the tunneling 
frequency, as is usually the case. The frequency shift 
u x = n g3 _ s (ve) is determined by the same average. 

The collisional master equation Q thus reduces to 
an equation in the 2d space spanned by the |V>f) = 
(|L) +e i7r ^|i?)) / V2- Denoting the tunneling frequency 
by uj z and identifying the energy eigenstates \ipo), 
with the eigenvectors of the Pauli matrix a z , it takes the 
form 

I 7 
dtp = 7T V°z6z + u x cr x ,p) + - (a x p& x - p) . (3) 
2,1 A 

This equation has the remarkable property of stabilizing 
the enantiomcr states \R) and \L), provided 7 3> w z , i.e., 



when the decoherence rate is much greater than the tun- 
neling rate. The chiral configuration states decay with 
the suppressed rate w^/7 <§C 7, w z in this case, while the 
superposition states |^) decay at least with rate 7 

The parity- sensitive dispersion interaction. — We now 
turn to the general van der Waals interaction between 
two polarizable particles to identify the relevant part that 
distinguishes a left-handed from a right-handed config- 
uration and thus gives rise to different scattering am- 
plitudes and f( R \ The total London dispersion 
interaction is conveniently expressed in terms of the 
frequency-dependent multipolc-polarizability tensors of 
the particles, most prominently those involving virtual 
electric dipole (ED) and electric quadrupole (EQ) tran- 
sitions Q . It is dominated by the standard van der Waals 
ED-ED/ED-ED interaction, which depends on the inter- 
particle distance r as r~ 6 and is determined by the ED- 
ED polarizability tensors a(iuj). This bulk interaction 
cannot distinguish different chiral configurations since 
the a (iuj) are parity-invariant. 

Among the chirally sensitive parts of the London 
dispersion interaction between a chiral and an achi- 
ral molecule, the most important contribution is due 
to the EQ-ED/ED-ED polarizability combination. It 
has a r _7 dependence and it combines the EQ-ED po- 
larizability tensor Aijk (iuj) of the chiral molecule 
with the ED-ED polarizability of the spherical projec- 
tile. One can safely neglect the contributions involving 
higher-order electric polarizability tensors due to their 
short-ranged nature; among the contributions involving 
magnetic susceptibilities, in particular due to magnetic 
dipole (MD) and magnetic quadrupole (MQ) transitions, 
the ED-MD/ED-MD vanishes if one molecule is achi- 
ral, while the ED-MD/ED-MQ and ED-MD/EQ-MD are 
suppressed by the square of the fine structure constant 
[ToL [llj]. Also the (parity- invariant) influence of a small 
permanent dipolc-momcnt of the chiral molecules can be 
usually disregarded, since it gets overshadowed by the 
van der Waals interaction. 

We note that the chirally sensitive EQ-ED/ED-ED in- 
teraction contribution is usually not accounted for in 
molecular scattering calculations. This is because its chi- 
ral disparity vanishes when averaged over all orientations 
of the chiral molecule, implying that the cross sections 
of left- and right-handed molecules are equal. However, 
such rotational averaging of the coupled-channel equa- 
tions, which drastically reduces the calculational effort, 
is not allowed in the present case, since the decoherence 
cross section ([1]) crucially depends on the phase differ- 
ences of the individual scattering amplitudes and 
. This phase information is lost at the level of stan- 
dard scattering cross sections, which indeed remain equal 
for left- and right-handed molecules. 

Solving the coupled- channel equations. — In order to 
evaluate the scattering amplitudes one must calculate the 
scattering matrix S by solving the associated coupled- 
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channel equations [12|]. Our system comprises a chiral 
asymmetric top molecule colliding with a noble gas atom, 
such that each channel is described by a set of rotational 
(pseudo) quantum numbers j, rnj, r, as well as by the 
orbital and the total angular momentum, £ and J. 

The resulting system of differential equations involves 
infinitely many closed channels and a finite number of 
asymptotically free, open channels, which scales with the 
total energy as E 3 / 2 . When truncating this infinite sys- 
tem, care is needed to ensure that at least all those closed 
channels arc retained which have an appreciable coupling 
to the initial state, even if their asymptotic energies lie 
deeply in the energetically forbidden domain. 

Since state-of-the-art numerical program packages for 
scattering calculations do not accommodate asymmet- 
ric molecules with their parity-changing collision dynam- 
ics, we developed a numerical method for solving the full 
quantum mechanical scattering problem (avoiding pre- 
mature rotational averaging). It is based on the log- 
derivative algorithm of Johnson and it allows to 
ensure the convergence with respect to the truncation 
of closed channels by adjusting threshold values both for 
the closed channel energies and for the internal angular 
momenta to be retained. We checked our program for 
the case of symmetric top molecules against the numeri- 
cal package Molscat [lj . 

Model for D2S2 ■ — The transition from tunneling to sta- 
bilization is best observed with a molecule of moder- 
ate tunneling frequency u) z . Motivated by recent pro- 
posals for enantiomer discrimination [TBI [l6| . we fo- 
cus on D2S2, one of the simplest chiral molecules, with 
lu z /2it = 176 Hz. Since ab initio calculations for the dy- 
namic tensors a (iui) and Aijk (iw) are not yet feasible for 
dihydrogcndisulfide (nor for any other chiral molecule) 
we extend the bond increment method (l7| for the calcu- 
lation of static molecular susceptibilities to the dynamic 
case. As described in [It}, static molecular susceptibil- 
ity tensors can be calculated as a sum over bond incre- 
ments associated to the constituent atoms a, which ac- 
count for the relative positions and orientation of the 
atoms and their mutual bonds. The bond increments are 
obtained from ab initio values of static molecular sus- 
ceptibilities of training compounds, which include the 
D2S2 molecule [l7j |. In a simple extension of this ap- 
proach, we replace the static values of the constituent 
atomic susceptibilities by a Drude model, a' a ) (iui) — 
aW(0)/„M and A\% (iuj) = A^ k (0) f a (u;), with 
f a (ui) = uj 2 / (lu 2 + uj 2 ), where ui a is the first excitation 
energy, ton = 0.375a.u and ujg = 0.252a.u. [Isj ]. The con- 
figuration of D2S2 is characterized by the bond lengths 
r ss =2.05A, r S D=1.34A, the inter-bond angle Z(DSS) = 
100.4° and the dihedral angle Z(DSSD) = 90.3°, and 
constants of inertia given in We take the back- 

ground gas to be ground state helium, whose (spherically 
symmetric) dipole-dipolc polarizability is accurately de- 




FIG. 1: Scattering cross section <T to t (upper panel, solid line) 
and decoherence cross section rjtot (lower panel, solid line) for 
the collision of a He atom off a ground state D2S2 molecule, as 
a function of the kinetic energy E/ks ^ 10K on a logarithmic 
scale. The dashed line corresponds to @, the dash-dotted 
lines gives the high-energy behavior ([5]). 

scribed by a sum of four Lorentzians [20| . 

Numerical results. — It takes substantial numerical ef- 
fort to observe the onset of collisional stabilization, since 
this regime is characterized by a large number of par- 
tial waves, while the standard semiclassical approaches 
for cross sections cannot be applied. Moreover, the num- 
ber of relevant scattering channels proliferates; at the ki- 
netic energy E/k-g = 300 K of the relative motion about 
1.2 x 10 4 coupled differential equations would have to be 
solved simultaneously, with as many different initial con- 
ditions, for extracting the S-matrix in the subspace of a 
single total angular momentum J. 

The upper panel of Fig. [T] presents the exact total 
scattering cross section a tot (for kinetic energies up to 
10K), starting from the rotational ground state of D2S2, 
while the lower panel shows the decoherence cross section 
7y to t = ^aVaa ■ One observes that a non-negligible effect 
of collisional stabilization can be expected already for en- 
ergies well below the threshold for the first channel that 
directly couples different parity subspaces (i^thres/^B = 
17.5 K). The decoherence cross section ij t ot tends to sat- 
urate at about 100a§, corresponding to Jftot/otot — 10% 
at 10K — a remarkably large value, given that the chiral- 
ity distinguishing interaction contributes only weakly to 

Ctot- 

High-energy approximation. — In order to understand 
this saturation and to assess the stabilizing effect of col- 
lisions at larger temperatures we now consider a high- 
energy approximation for n to t and Utot- Since the Born 
approximation renders the scattering amplitude a lin- 
ear functional of the potential, it follows immediately 
from (fTJ) that the asymptotic high-energy behavior of 
rjtot is the cross section corresponding to Ay (r) = 
Vl ( r ) — Vr ( r )- We start from the general dependence of 
the cross section on the initial state elastic S-matrix ele- 
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ment, a tot = 2tt £j (2 J + 1) (1 - Rc((* |^( J ) |* )))A 2 
[l2| . apply the exponential Born approximation, and take 
the high-J asymptotics. For a homogeneous, spherical 
potential V (r) = C n r~ n , this yields the standard result 

(7 = p n {Cf^/E) 1 ^ n 1 , with p n a numerical factor de- 
pending weakly on n > 2 12 1 . 

Disregarding the molecular core, we can thus approxi- 
mate Ctot by the dominant, spherically symmetric n = 6 
contribution to the van der Waals ED-ED /ED-ED inter- 
action, i.e., (Ttot — CTjot with Cq = 11.7a.u. Also AT^ (r) 
is homogeneous in r, with n = 7, but there is no di- 
rect Born contribution to the decoherence cross section 
?7tot since (^ol^F (t) l^o) vanishes. Therefore, we in- 
clude all off-diagonal elements that couple the ground 
state in the exponential Born expression. This yields 

ReC^ol^l^o)) * cos(^E/ #0 K^olAyl*/)! 2 ), with 
^ / the free channel wave functions in the J subspace. 

The lowest channel V^ithat couples to 4"o via AF opens 
at 17.5K, with j = 3. Expressing the corresponding 
coupling as |(* |AF|*i)| = h 2 P 5 / (2m»r 7 ) , with (3 a 
convenient parameterization and m* the reduced mass, 
and replacing the sum by an integral this yields 
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with i? = (tt) 2 /2m,. This approximation is displayed 
as the dashed line in the lower panel of Fig. [TJ The com- 
parison with the exact results indicates that (O starts 
to apply already at numerically accessible energies. An 
asymptotic expansion yields 



Vtot — ci 
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fc 7/6 



(5) 



with ci = 3.66 and C2 = 14.4 (see dash-dotted line in 
Fig. [J). 

Discussion. — Our numerical and analytical analysis 
suggests that the ratio rj tot /a tot varies only weakly at 
larger energies. A typical value at 300 K of r} tot /cr tot ~ 
25% means that environmental stabilization sets in once 
the rate of collisions exceeds four times the tunneling 
frequency. Noting that D2S2 tunnels with lo z /2h — 
176 Hz, we find from ((5]) that the critical pressure 
and temperature of a helium atmosphere must satisfy 
(p/mbar) (T/K) _2/3 > 3.0 x 10~ 7 , i.e., p ^ 1.6 x 
10" 5 mbar for T = 300 K. 

This prediction could be tested in an experiment that 
applies laser based coherent control techniques in a Stcrn- 



Gerlach-type setup for separating a molecular beam into 
left- and right-handed daughter beams 15|, 1(|. Passing 
one of them through a gas cell, one may analyze the enan- 
tiomeric purity with another Stern-Gerlach stage. Since 
thermal racemization can be controlled by cooling (given 
the 2300 Kfcs barrier height for chiral flipping in D2S2), 
it is thus possible to directly observe the expected en- 
vironmental stabilization of the chiral configuration as a 
function of the gas pressure. 
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